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Abstract - Cosmological models in Lyra's geometry are constructed and investigated 
with the assumption of a minimal interaction of matter with the displacement vector 
field and the dynamical A term. Exact solutions of the model equations are obtained 
for the different equations of state of the matter, that fills the universe, and for the 
certain assumptions on the decaying law for A. 
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! 1 Introduction 



The study of the cosmological models with a varying A -term has been essentially 
activated by the efforts of several researchers from as long as a cosmological constant problem 
gained a real confirmation of decaying A(t). There are a large number of observations for 
the determination of Einstein's cosmological constant, A, or some matter ingredient of the 
Universe which is slowly changing with time and, thus, acting like A. Moreover, recent 
observations of supernovae of type la [TJ [5] , Cosmic Microwave Background Radiation [3J SJ , 
Baryon Acoustic Oscillations in galaxy surveys [5j[6] etc., are of evidences in favor of a non- 
zero cosmological 'constant' with relative energy density J1a = A/3-ffg « 0.6 — 0.7. This 
\ value could not remain constant during long-term observation. The recent studies of nonlocal 

effects, wormholes, inflationary mechanisms of cosmological perturbations are the evidence 
in favor of decreasing with time of effective cosmological term. Cosmological models with a 
■ time- varying A(t) with different decay laws were proposed by several researchers during the 

C last two decades (see, e.g. [7] -[H])- 

However, the formal introduction of the dynamical A - term into Einstein's equation of 
General Relativity (GR) leads to the violation of the energy conservation law of matter. This 
follows directly from the Einstein's field equation due to the Bianchi identity for the Riemann 
curvature tensor, if the cosmological term is constant or zero. Previously, several authors 
have investigated cosmological models in Lyra's geometry |15j . which is a generalization of 
Riemannian geometry by introducing a gauge function which removes the non-integrability of 
the length of the parallel transfer characteristic of Weyl's theory (see, eg, [16] and references 
therein) . It was noted that cosmology based on Lyra's manifold with constant gauge vector 
is similar to the C- field theory of Hoyle-Narlikar [IT, or contains a vacuum field, which 
together with the gauge vector field can be considered as a cosmological term. A scalar-field 
cosmological model in Lyra's geometry is investigated in |18j . where it was also noted that 
the dynamical displacement field, which is free from interaction with matter and does not 
violate the energy conservation law of matter, can only serve as a stiff perfect fluid. Thus, we 
can assume that the simultaneous consideration of the dynamical A - term and displacement 
field is able not only to prevent a violation of the energy conservation law but also to enrich 
the theory. In our view, this approach could lead to the significant variations in evolution 
of the standard cosmological model. 

The purpose of this paper is to construct a cosmological model in Lyra's geometry, pro- 
vided to the minimal interaction of matter with the displacement vector field and dynamical 
A - term. We find some exact solutions of the model equations for the different states 
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of matter that fills the universe, and certain assumptions regarding the evolution of the 
cosmological term. In addition, we analyze the possibility of the accelerated cosmological 
expansion on the basis of solutions obtained here. 



2 A Model with Time Evolving Cosmological Term in 
Lyra's Geometry 

The Einstein's field equations with a cosmological A - term in Lyra's geometry, as proposed 
in |19j in normal gauge, may be written as 

1 3 3 

Rik - -xgikR - A^ife + 2 ~ -^Sikfj*' <t>j = To,, (l) 

where fa is the displacement vector, the gravitational constant is 8ttG = 1, and other 
symbols have their usual meanings in the Riemannian geometry. The energy-momentum 
tensor (EMT) of matter Tik can be derived in a usual manner from the Lagrangian of 
matter. Considering the matter as some effective perfect fluid, we can write: 

Tik = (Pm + Pm)UiUk ~ Pm 9ik, (2) 

where U{ = (1, 0, 0, 0) is the 4-velocity of the co-moving observer, satisfying UiU 1 — 1. Then 
let us represent fa as a time-like vector field of displacement, 

-^=AO,0,oY (3) 

where ft = f3(t) is a function of time alone, and the factor 2j\[Z is introduced to simplify 
the subsequent equations. The metric of a Friedmann-Robertson- Walker (FRW) space-time 
can be written as follows: ds 2 = dt 2 — a 2 (t)(dr 2 + £ 2 (r)d£l 2 ), where a(t) is a scale factor 
of the Universe, £ (r) = sin r, r, sinh r in accordance with the sign of the spatial curvature 
k = +1, 0, — 1. Given this metric and Eqs. d2j) , {J3J) , the field equation ([!]) can be reduced to 
the following set of equations: 

3fc ■ k 

3H 2 + ^-/3 2 =p m +A, 2H + 3H 2 + ^ + /3 2 =- Pm +A, (4) 

where H — a/ a is the Hubble parameter, and an overdot denotes differentiation with respect 
to time t. 

As a consequence of Eq. the continuity equation for the effective matter can be 
written as: 

p m + A + 2(3 p + 3H [p m + p m + 2/3 2 ] = 0. (5) 

One of the most important quantity to describe the features of dark energy models is the 
equation of state (EoS) It is known that the most important quantity to describe the possi- 
bility of an accelerated mode of expansion is the so-called deceleration parameter q, defined 
as 

To proceed further in studying of our model, it is necessary to determine the type of 
dependence of the cosmological term (or the displacement vector field) on time, or the type of 
interaction between the geometric fields and matter. For simplicity, we consider a spatially 
flat FRW cosmology with k — 0. We can write the basic equations of the model ((4]) as 
follows: 

3H 2 = Pm + A + p 2 , 2H = -( Pm + Pm + 2f3 2 ), (7) 

and the continuity equation ([5]) is in the same form. The latter follows from the set of 
equations ([Tj) . 
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The main assumption on our model consists of the minimal coupling of the matter with 
the displacement vector field and the cosmological term on Lyra's manifold. That means 
that there is no the direct interaction between them, that is this interaction is realized only 
through the gravitational field. It allows us to avoid the violation of the energy-momentum 
conservation of matter in the framework of minimal demands on the behavior of displacement 
vector. Indeed, due to the covariant equation of the energy-momentum conservation T* fc = 
and the identity G\. k = for the Einstein tensor, the field equation (TJJ leads to the following 
equation for the cosmological term and the displacement field: 



|(M fc -^&) ifc =0. (8) 



In the absence of cosmological term (or in the case of its constancy), when A ; j = 0, Eq. 
(JU) leads to the interpretation of the displacement field as an analog of the so called stiff 
fluid with the EoS w s f =1. The assumption of the non- vanishing (and varying with time) 
cosmological term can significantly extend the capabilities of such a model in describing of 
the real processes taking place in the Universe. 

Here, the minimal coupling means that the energy conservation law for matter is valid 
regardless of the presence of A(i) and f3(t): 

p m + 3H(p m + p m ) = 0. (9) 

Thus, the remaining part of Eq. ([5]), i.e. Eq. ©, has the form: 

A + 2/3/3 + 6H/3 2 = 0. (10) 

The set of equations ([7]), ©, (fTU)) . which determines the dynamics of our model, must 
be supplemented by some conditions. In our view, the most realistic conditions are: first, 
the EoS of matter p m — w m p m and, second, the law of evolution of the cosmological term 
A(t) which should correspond to the observational data. Of course, this case is not the only 
one in searching of exact solutions in the framework of such a model. Nevertheless, in the 
present paper we study our model just with this assumption, as we have no any possible 
dependence j3{t), which could be proved by the direct observations. 

We assume that the matter content of the Universe can be considered as a barotropic 
perfect fluid with a constant equation of state (EoS): — 1 < w m < 1. Then Eq. (|9]) can be 
easily integrated which yields 

p m ^ Po a-^ 1 + w ^, (11) 

where po is a constant of integration. Due to (jTTJ) , the equations (O can be written as 
follows: 

■AH 2 = p a _3(1 + Wm) +A + /3 2 , (12) 

2H = -(l + w m )p„a- 3 ( 1 + w ^ -2/3 2 , (13) 



and Eq. (flQ|). as a differential consequence of Eqs. (fT2j) . ([T3|) . is not changed, that can be 
easily verified. 

Thus, the problem is reduced to solving the set of equations (|12j) . (fl3|) either for H(t) 
and /3(i) with some given function Kit), or for Hit) and A(i) with a given function Pit). 
Since no information about the possible dependence of /3(t) is available, we consider the 
phenomenological models with some given functions A(£) that have an observational ground 

and are widely discussed in the literature (see, e.g. [H]), namely: it) A(t) = — , and 

t 2 

iii) A{t) — —Hit). Besides, two cases, w m — —1 and w m ^ —1, should be considered 
separately. 
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Figure 1: Plots of a(£) and Hit) vs time 
for the w m = — 1 model in the case A e f / = 
a/t 2 . Here to = 1. 



Figure 2: The plot of g(£) vs time for the 
A e // = a/t 2 model. 



3 Models with the EoS of quasi-vacuum: w r 
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In this case, Eq. (jlip leads to the usual property of the quasi-vacuum state: p m = —p m = 
Pv = constant. Then, Eqs. (1121) and (|13|) can be rewritten as 



H + 3H 2 = p v + A, (3 2 



H. 



(14) 



It is easy to see that the constant energy density of the quasi-vacuum po can be added to 
the cosmological term. By introducing notation A e ff = p v + A, we can obtain from Eq. 
(fl4)) : H + 3H 2 — A e fj. We consider two cases depending on A e //(£), mentioned above. 

As a result of solving Eq. (fT4f with A e ff — a/t 2 for the case (i), we obtain the following 
expressions: 



H(t) 
a(t) 
(3 2 (t) 



— { 1 + n a tanh 
6t I 



■y ln(£/£ )J } , n a = y/1 + 12a, 
£ > 
£o ' 



1 



12£ 2 



1 + n„ tanh 



ii 



rHt/to] 



(15) 
(16) 
(17) 



The plots of a(t) and H(t) for this solution with two different values of the coupling constant 
a are shown in Fig. 1. 

It is easy to find that the deceleration parameter © is equal to 



q(t) =2-3- 



1 



tanh 



ln(£/£ )]} 



(18) 



due to Eq. (fT5|) . The plot of q(t) for various values of the coupling constant a is shown 
in Fig. 2. One can see that for each value of a there exists a finite positive value t cr , at 
which the deceleration parameter is the negative infinity. This t cr can be found from the 
equality of the denominator in Eq. (fT5|) to zero. Besides, from Eq. (fT7| it follows that 
the displacement held becomes real only since a certain instant £j, which can be found from 
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Figure 3: Plots of a(t) and H(t) for the 
A e ff = (a/t)H model with ti = l. 



Figure 4: The plot of q vs time in the 
A e // = (a/t)H model for some values of 
a and U = 1. 



equation j3 = 0. Thus, we have 



1 



n a + 1 



Iq, - 1 + a/?^ 



1 



(19) 



One can see that ij = for n a = 1, i.e. for a = 0. The comparison of formulas in Eq. (|T9|) 
leads to the conclusion that, the inequality U > t cr is valid, except the trivial case a = 0. 
Moreover, for < a < 1/12 we have < to, and we have always U > to for a > 1/12. 
These estimates suggest that the nature of this model becomes realistic when t geqti, that is 
when there are no any problems with the displacement field to be real, and with the infinite 
acceleration of expansion. 

For the case (ii), we get from Eq. (fT4| the following equation for the Hubble parameter: 



H 



ZH Z = -H, 



whose solution can be easily found as follows: 

a(l + a) 



H(t) 



{t/U 



3fe 



1 + a(t/t i ) Q+1 

where t% is a constant of integration. Then integrating Eq. 
factor in the following form: 

r i V 3 

a(t) =a l + a(t/tj) 



(20) 
(21) 

TJ for a(t), we find the scale 
(22) 



where clq is a constant. Taking into account the second equation in (|14p and Eq. (|21[) . the 
following expression for the displacement field can be found: 



(t/U) 



a+l 



1 



(23) 



[l + a(t/t,) Q+1 ] 2 

From this function (3 2 (t), it can be seen that (3 2 changes its sign at the instant U > 0. In 
view of the Hubble parameter (I21[) . we can obtain by simple mathematical manipulation 
that the deceleration parameter © is as follows 

3 



lit) 



1 



1 



1 



(24) 
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The plots of the scale factor a(t) and the Hubble parameter H(t) in this case is shown in 
Fig. 3. The plot of q(t) versus time is shown in Fig. 4. 



4 Models with the EoS of matter w m ^ — 1 

Multiplying Eq. (fl"2|) by (1 + w m ) 7^ and adding the result to Eq. (fT3|) . we obtain the 
first equation of the set of independent equations in the following form: 

2H + 3(1 + w m )H 2 = (1 + w m )A - (1 - w m )(3 2 . (25) 

We can take Eq. (|T0|) as the second independent equation. 

For this case, we replace the assumptions on A(t) considered above by assuming that the 
cosmological term is proportional to the invariant of the displacement vector field. Hence, 
in this section we are supporting the idea expressed earlier about the possible role of the 
displacement field as an effective A - term in the Einstein equation: A = 7/3 2 , where 7 is an 
arbitrary positive constant. Then Eq. (|10p can be integrated by introducing a new variable, 



In 



dx 

a{t)/a{t m ) H{t) = — , (26) 



which coincides with the number of e-folds and means expansion of the universe in e x times 
at present versus the scale factor o, = a(ti n ) at the initial time U n . As a result, we have 
that 

6 6 



A(x) = A e 1+7 ,/3 2 ( x ) = ^ e 1 + 7 , (27) 

7 

where Ao is a constant of integration. Note that it follows from Eq. (|2"T|) due to Eq. ([2"B)) 
that the following law is realized: A(t) oc a(t)~ b , where b > 0. This decay law for the 
variation of the cosmological term was investigated earlier, for example, in [7], [20] . 
Substituting Eq. (j2"7| into Eq. (|2"5)l . we obtain the following equation: 

6 

dH 2 ~1 x 

" +3(l + w m )H 2 = A [(l + w m )-(l-w m )j- 1 ]e 1 + 7 . (28) 



-.r 



dx 

The general solution of this equation can be written as 



(dx\ 2 _ A —bx 1 r<„—Bx 



where 



H '^Tt)=— b + (29) 



A = Ao[(l+w m )-(l-w m )7 _1 ] = 2A 



6 - b ' 



B = 3(1 + w m ), b=^—; B, b G (0; 6), (30) 

1 + 7 

and C is a constant of integration. In this case, as seen from Eq. (|2T))) , we assume that 
B 7^ b, i.e. 7 7^ 70 = (1 — w m )/(l + w m ). Otherwise, i.e. for 7 = 70 +> A = 0, the general 
solution of Eq. P5| can be written as: 

H 2 = (IE.) 2 =Ce -3(l + w m )x (31) 
V at J 

where C > 0. Taking into account the definition (|2"B)) . it is easy to verify that the solution 
of Eq. (|31l) reproduces the standard solution for FRW model: a(t) oc t 2 / 3 '^ 1+Wm \ 

In the case of 7 7^ 70 , the general solution of Eq. (|29[) can be written only in an implicit 
form as 

a(t)/a,i 

t T - d * =±J^(l + 1 -)t + C 1 , (32) 

Vc + Z B - b V 3 V 1) 
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W m = - 1/2 — ■ — w m = - 1/3 w,„ = 1/3 w,„ = O 



Figure 5: Evolution of the deceleration parameter in the A = 7/3 2 model for the different 
EoS of matter w m =/= —I with 7 = 5 and Co = 1. 



where Co = C(B — b)/A is a dimensionless constant, C\ is an integration constant, and the 
equation A/(B — b) = A (l + 7)/37, following from (|30|) . is used. The integral in Eq. (|32j) 
can be found for some specific values of B and b, defined by the expressions (I30[) through 
the EoS of matter w m and the coupling constant 7. So, say for w m = =>■ £> = 3 and 
7 = 7 =*> = 3/4 in (021), we have: 



a(t) 



= — 1 + - cosh 



8A 
21 



where we put C\ = 0. It follows that q(0) — 1/2 > at the initial time, and then at 
some moment to 5 defined by equation q(to) = 0, the decelerated expansion followed by the 
accelerated one: q(t > t ) < 0. 

However, the model can be studied for the presence of acceleration in the general case 
with the help of definition (|6]) and Eqs. (|25[) . (f27|) and (|29|) . As a result, we can obtain the 
following expression: 

B (B- b)e( B -V x 
2 "> l"c + e( B -V x 



q(x) 



-1 



(33) 



where B ^ b 7 7^ 70. As follows from this expression, the deceleration parameter evolves 
from qi = q(—oo) = —1 + B/2 at the beginning of expansion up to qf = g(+oo) = —1 + 6/2 
to date, provided to B > b, but - from qf to qi for B < b. Note that the denominator in the 
formula (|33l) is non-negative, as it follows from the right hand side of Eq. (|29[) rewritten in 



the form H 2 = (A /3)(l + 7- 1 ) e 



Bx 



Co 



D (B-b)x 



] > 0. Since the latter must be true for 



all values of x € (—00, +00), it also follows that Co > 0. The evolution of q(x) according to 
Eq. (|33[) for several values of w m and Co = 1 is shown in Fig. 5. 



5 Conclution 

Thus, we have studied the cosmological models in Lyra's geometry, supposing the so- 
called minimal interaction of matter both with the displacement vector field and with the 
dynamical A - term. Besides general study, we have given some examples of exact solution 
for the model under consideration. Exact solutions of the dynamical equations for our model 
in the cases w m = — 1 and w m =^ — 1, and under various assumptions about the evolution 
of the cosmological term are obtained. Surely, these solutions can not cover all possible 
applications of the present research. 

Nevertheless, an interesting feature of these models with a quasi-vacuum EoS. As can be 
observed in Fig. 2 and Fig. 4, these models begin to expand with a super-acceleration, and 
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then asymptotically approach a state of the constant acceleration or the non-accelerated 
state of expansion with q > 0, depending on the value of the coupling constant 7. It could 
be assumed that these models are relevant to the cosmological inflation. 

Cosmological models, built in the last section with the assumption of proportionality 
between the cosmological term and the displacement field, demonstrate a tendency to shift 
from the slow expansion to the accelerated expansion. That can be seen in Fig. 5. This be- 
havior of these models suggests its relation to the phenomenon of the late-time acceleration, 
which is reliably confirmed by the observational data. 

The interesting idea, which could be realized, consists of the combination of two cases 
considered above. Indeed, the behavior of the deceleration parameter in that case might 
reveal the inflationary acceleration in the very beginning of expansion as well as the late- 
time acceleration. Moreover, it is easy to obtain the exact solution for our model under 
the assumption of some different functions A(i) widely discussed in the literature, namely 
A cx H 2 (t), A cx aH + -yH 2 etc. Further details and consequences of the model considered 
here are in progress. 
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